Cell Reports

Wildfire smoke is an escalating global public health threat exposing millions of people, including children, to hazardous air pollution each year. Although wildfire smoke toxicants have been linked to a range of adverse health outcomes, including immune dysregulation, the long-term consequences of real-world pediatric wildfire smoke exposure on health and development remain largely unknown. To investigate the persistent effects of early-life exposure on immune health, here we leveraged a cohort of rhesus macaques that experienced nine consecutive days of hazardous wildfire smoke exposure in infancy during the 2020 Oregon Labor Day wildfires. By integrating ex vivo immune stimulations, multiplex cytokine profiling, single-cell transcriptomics, and genome-wide DNA methylation profiling, we identified persistent immunological consequences across molecular and functional levels. We found that a single severe postnatal exposure, in the first three months of life, was associated with persistent change in the innate immune response, including reduced pro-inflammatory cytokine response to a bacterial endotoxin, with subtle but consistent transcriptional changes in myeloid cells, particularly among males. Wildfire smoke exposure was also associated with changes in proportion of B and T/NK cells, and within the T/NK cell compartment, exposed animals exhibited an expansion of cytotoxic cells. Consistent with this, CD8+ T cells displayed extensive transcriptional remodeling and shifted toward more differentiated effector states, with the greatest differentiation observed in animals exposed at the youngest ages. Genome-wide DNA methylation profiling identified smoke-associated methylation changes consistent with acceleration of epigenetic aging, as well as persistent epigenetic alterations impacting genes involved in oxidative stress responses, innate immunity, T cell differentiation, and hematopoiesis. These findings demonstrate that a single severe wildfire smoke exposure during a critical developmental window is associated with extensive immune and epigenetic remodeling that persist years after exposure, providing new insight into the long-term biological consequences of early-life wildfire smoke exposure.

High-grade serous ovarian carcinoma (HGSC) responds poorly to immune checkpoint blockade, partly due to a macrophage-dominated immunosuppressive microenvironment. We integrated single-cell spatial proteomics and spatial transcriptomics across 50 HGSC tumors and applied SPACEstat to resolve higher-order immune communities and their transcriptional programs. We identified six immune community types, with macrophage-dominated Myelonets representing the predominant spatial pattern of immune organisation. In chemotherapy-exposed tumors, Myelonets showed coordinated lipid metabolism-immunosuppression and inflammation-MHC-II macrophage transcriptional programs, with SPP1, C1Q, VEGF, MMPs, and CCL18 linked to immunosuppressive states and fibroblasts emerging as key mediators of macrophage communication. Chemotherapy contracted large Myelonets while increasing CD8+ T-cell organization into Lymphonets. Persistent macrophage dominance within Myelonets was associated with adverse outcomes among patients who achieved a complete response to treatment. Together, we identify Myelonets as clinically relevant, multicellular immunoregulatory niches sustained by spatiotemporally coordinated macrophage programs and stromal crosstalk.

To sustain blood formation, hematopoietic stem and progenitor cells (HSPCs) coordinate a multitude of cell biological processes, from cell cycle control and stress responses to lineage priming. While many genetic regulators of high-level HSPC function have been identified, how HSPCs coordinate more basal cell biological programs, and how such programs relate to stem cell function, remains incompletely understood. Here we use Perturb-seq to profile the transcriptional consequences of targeting 520 genes by CRISPRi in primary mouse HSPC cultures. We developed an analytical strategy to separate perturbation-induced changes in cell-state abundance and clonal heterogeneity from cell-state-local transcriptional effects. From these local perturbation signatures, we identified 19 gene regulatory programs (GRPs) that are defined by co-regulation in response to genetic perturbation, in contrast to co-expression or human curation, and align well with cell biological processes. By decomposing gene expression data from functional and clinical studies into program activity, we show that GRP activities associate with, and predict, phenotypes such as clonal output after transplantation, as well as survival and drug response in retrospective acute myeloid leukemia (AML) cohorts. Together, our study establishes perturbation-derived co-regulation programs as an interpretable framework for linking genetic regulators, cell-biological processes and stem-cell-associated phenotypes.

Integrating multisensory and behavioral information is essential for sensory perception. In the auditory system, multisensory and behavioral influences emerge early in subcortical structures. Descending projections from non-auditory cortical areas are well positioned to convey such signals, yet how they shape subcortical auditory processing remains poorly understood. Here, we investigated corticocollicular projections from the primary somatosensory (S1) and motor (M1) cortices to the inferior colliculus (IC), a principal integration center in the auditory midbrain. We found that trunk- and limb-related regions of S1 and M1 form prominent monosynaptic projections to the IC, and that optogenetic activation of these projections robustly drives IC activity. Notably, a substantial population of cortical-responsive neurons did not respond to sound. In sound-responsive neurons, concurrent cortical stimulation enhanced sound-evoked responses, whereas cortical activation preceding sound onset suppressed them. Furthermore, both cortical-responsive IC neurons and deep-layer S1 and M1 neurons exhibited locomotion-related modulation and anticipatory activity prior to movement onset, suggesting that these descending pathways convey movement-related signals to the IC. Together, our findings identify a descending sensorimotor circuit that integrates body- and movement-related information with auditory processing in the auditory midbrain.

Dissecting the complexity of antibody responses in orthopoxvirus (OPXV) infected individuals is essential for elucidating protective mechanisms and identifying candidate protective immunogens. Here, we profiled the acute humoral response in 51 mpox cases, showing distinct IgG trajectories among multiple antigens alongside the rise of plasma neutralizing activities to plateau within 6 weeks after symptom onset. Utilizing a single-cell transcriptomic and BCR sequencing based antigen-agnostic mAb isolation workflow, we further generated monoclonal antibodies (mAbs) from 254 expanded peripheral B cell clones of 3 patients. We discerned 97 specific mAbs recognizing at least 12 different OPXV proteins via integrated screening approaches, which comprised neutralizing antibodies binding unconventional viral targets and antibodies exhibiting extraordinary in vitro and in vivo anti-OPXV effects. The number of OPXV-specific mAbs recovered per donor reflected the percentage of expanded clones among circulating B cells. More interestingly, we demonstrated that the inferred unmutated common ancestors (UCAs) of neutralizing antibody clones did not necessarily react with OPXV, implying that OPXV neutralizing antibodies might frequently originate from B cells previously activated by unknown antigens. Our work establishes an efficient workflow for antigen-agnostic isolation of pathogen specific mAbs and reveals previously unclarified features of antibody responses induced by acute MPXV infection.

Trisomy 21 or Down syndrome (DS) affects multi-organ systems across the lifespan. The presence of an extra chromosome, along with genome dosage imbalance due to triplicated genes, contributes to the DS phenotypes. Of the DS mouse models, few are aneuploid with a freely segregating extra chromosome. We previously showed that the aneuploid Ts65Dn mice exhibit metabolic deficits consistent with the metabolic profile of DS. However, the genotype-phenotype relationships in Ts65Dn mice are complicated by the presence of triplicated genes unrelated to human chromosome 21 (Hsa21). To address this issue, we leveraged a refined model, Ts66Yah, where the extra triplicated genes in Ts65Dn have been removed. Deep phenotyping and multi-omics analyses showed that Ts66Yah mice develop pronounced and widespread metabolic disturbances. Despite sexual dimorphism in weight gain, body temperature, lipid and lipoprotein profiles, hepatic injury and adipose fibrosis, both male and female Ts66Yah mice share a common phenotype of pronounced glucose intolerance and insulin resistance, reduced mitochondrial respiratory capacity in visceral fat, altered serum inflammatory cytokine profile, and dysregulated serum and liver metabolomes. Pan-tissue transcriptomes also reveal signatures of immune activation, disrupted metabolic processes and cellular respiration, altered cytokine signaling, enhanced oxidative stress, and extracellular matrix remodeling. These combined changes across tissues disrupt metabolic homeostasis more severely in Ts66Yah than in Ts65Dn mice. Several phenotypes, including glucose intolerance, insulin resistance, tissue fibrosis, and oxidative stress were further exacerbated by an obesogenic diet. This foundational data establishes Ts66Yah as a valuable reference model for the mechanistic and comparative study of metabolic dysfunction in DS.

Increased activity of leucine-rich repeat kinase 2 (LRRK2) confers Parkinson's disease risk. LRRK2 dynamically localizes to lysosomal membranes in response to various stresses, yet the mechanisms by which distinct lysosomal perturbations are communicated to LRRK2 remain unclear. Here, we show that inhibition of the lysosomal lipid kinase PIKfyve promotes LRRK2 recruitment and signaling through a pathway that requires the lysosomal chloride/proton antiporter ClC-7. ClC-7 in turn controls the accumulation of multiple Rab GTPases on lysosomes. LRRK2 signaling under these conditions requires its established Rab-binding surfaces, with Rab12 contributing significantly to this response. This pathway operates independently of CASM. In contrast, lysosomal stresses that induce CASM require both Rab-binding sites on LRRK2 and GABARAP for robust LRRK2 signaling. These findings identify ClC-7-dependent lysosomal remodeling and Rab accumulation as key features linking PIKfyve inhibition to LRRK2 signaling and reveal that distinct lysosomal stresses engage different combinations of Rab and GABARAP inputs to activate LRRK2.

Nonarteritic anterior ischemic optic neuropathy (NAION) is the leading cause of acute optic nerve related vision loss in older adults, yet no disease modifying therapy exists. Although ischemia is a defining feature of NAION, prior therapeutic efforts targeting vascular insufficiency or nonspecific oxidative stress have failed to prevent irreversible retinal ganglion cell (RGC) degeneration, underscoring an unresolved mechanistic gap between ischemic insult and permanent axonal failure. In this endeavour, we identify lipid peroxidation as an important driver of neurodegeneration in NAION. Analyses of human NAION retina, together with a rigorously validated mouse model, demonstrated a remarkable activation of phospholipid peroxidation within the retina following ischemic injury. RGC-specific overexpression of glutathione peroxidase 4 (GPX4), the only known enzyme capable of directly detoxifying phospholipid hydroperoxides within biological membranes, confers striking protection of RGC survival, axonal integrity, and visual function. We further demonstrate that mitochondrial-targeted GPX4 provides superior protection, suggesting mitochondria as a critical locus of lipid peroxidation-driven vulnerability in NAION. Leveraging real-time multiparametric in vivo imaging to directly interrogate axonal metabolism and function, we demonstrate that RGC-specific GPX4 overexpression robustly restores axonal and retinal mitochondrial abundance, improves ATP bioenergetics, and suppresses superoxide stress following optic nerve ischemia. Mitochondria-targeted GPX4 expression further restores axonal transport and retinofugal projections to central visual targets, thereby stabilizing visual pathway connectivity. Notably, these neuroprotective effects are recapitulated by Ebselen, a clinically tested GPX mimetic, identifying lipid peroxide detoxification as a translatable and imaging-validated therapeutic strategy. Collectively, this work establishes ischemia-induced lipid peroxidation as an essential driver of neurodegeneration in NAION and identifies GPX4 as a key molecular determinant of retinal ganglion cell resilience.

Post-acute infection syndromes (PAIS) follow viral syndromes including post-acute sequelae of COVID19 (PASC) which complicates 10-25% of SARS-CoV-2 infections. These syndromes lack precise explanatory mechanisms. We studied 173 human saliva proteomes during respiratory viral syndromes, seeking associations between 44 clinically-relevant protein expression patterns and subsequent sequelae counts. Exploratory models adjusted by extensive clinical annotations found interactions between 23 acutely-responsive proteins and SARS-CoV-2 infection that inversely predicted subsequent neurocognitive sequelae. An overlapping 19 acutely-responsive proteins during any acute respiratory viral syndrome inversely predicted general fatigue-related sequelae. Altogether, 29 proteins, derived from interferon stimulated genes (ISG), were uniformly beneficial, including 13 predictive of both neurocognitive and general sequelae. The proteins suggested both shared early pathobiology and virus-specific protective responses that shaped resolution of acute disease and different PAIS. Acutely elevated protective ISG proteins associated with reduced post-viral symptoms identify investigational starting points for novel mechanisms, diagnostics and therapeutics for PASC and PAIS.

The QxxR motif is evolutionarily conserved within DEAD-box RNA helicases, including Drosophila Me31B and human DDX6, which post-transcriptionally regulate gene expression during animal development. A pathogenic H372R substitution (QxHR to QxRR) in the QxxR motif of human DDX6 has been associated with various developmental defects, but how this motif contributes to DDX6-family protein function remains unclear. Here, we used Drosophila Me31B as an in vivo model to investigate the QxxR motifs developmental role. We generated a Drosophila strain carrying the corresponding H333R missense mutation in Me31B and characterized its effects on female fertility, embryonic viability, germline development, and Me31B-associated molecular pathways. The me31BH333R mutation reduced female fertility in a gene dose-dependent manner, with homozygous mutant females being sterile. Embryos from the mutant females also exhibited primordial germ cell defects. Despite these developmental phenotypes, the me31BH333R mutation did not significantly alter Me31B protein abundance, global ovarian transcriptome or proteome profiles, or representative germ plasm mRNA and protein localization. In contrast, bait-normalized IP-MS analysis revealed altered enrichment of selected Me31B-associated proteins, including increased association of known Me31B interactors Trailer hitch (Tral) and Ypsilon Schachtel (Yps). These findings establish Me31BH333R as an in vivo model for investigating the conserved QxxR motif and suggest that disruption of this motif compromises development not through broad changes in gene expression, but potentially through altered composition or regulation of Me31B-containing ribonucleoprotein complexes.

Microorganisms across all domains of life engage in molecular conflict, deploying toxins to inhibit competitors or respond to biological threats. Among these, ribonuclease toxins are particularly widespread and diverse. A substantial fraction is associated with the BECR fold, a compact /{beta} architecture that supports RNase activity despite extensive divergence. Although several canonical members are well characterized, many BECR-fold proteins remain difficult to identify because of low sequence similarity, variation in catalytic residues, and structural elaborations that obscure evolutionary relationships. The growing availability of high-confidence protein structure predictions provides an opportunity to reassess this deeply divergent protein landscape. Here, we integrate iterative profile-HMM searches, profile-similarity networks, structural analyses, active-site mapping, and genomic context to examine BECR proteins across the tree of life. Our analysis resolves an expanded BECR-fold landscape comprising canonical BECR and BECR-like superfamilies, refines the organization of canonical BECR proteins and identifies previously unrecognized families. We further validate BECR-Tox2 as a toxin neutralized by a cognate immunity protein and show that its homologs occur in both Menshen-like anti-phage systems and polymorphic toxin loci. Together, these findings expand and clarify the BECR-fold landscape and provide a framework for identifying and interpreting highly divergent proteins of this fold.

Vaccine responses vary across individuals due to differences in ageing and health status. Using transcriptomic profiling, we analyzed early gene expression profiles after influenza (QIV) followed by pneumococcal (PCV13) vaccination in 148 participants spanning young, middle-aged, and older adults. The two vaccines induced distinct immune signatures: QIV elicited innate and interferon immune activation, while PCV13 triggered inflammation-based responses. Older adults showed weaker but similar transcriptomic profiles compared to young adults. Among older adults, frailty, in addition to age, was strongly associated with reduced innate responses. In addition, we identified associations between early-stage transcriptomic profiles and later-stage antibody responses for QIV; however, no such associations were observed for PCV13. Importantly, observed group differences arose not from altered immune modules but from differences in the magnitude of gene expression, paving the way for immune-boosting interventions to enhance early gene expression in at-risk populations.

Clonal hematopoiesis (CH) increases with age, but whether different somatic clones represent an ageing phenotype or exert distinct systemic effects is unclear. In 450,587 UK Biobank participants, including 46,324 with plasma proteomics, we compared clonal hematopoiesis of indeterminate potential (CHIP) and mosaic loss of chromosome Y (mLOY) or X (mLOX) across biological ageing, incident disease, and circulating proteins. Despite shared age dependence, these alterations showed distinct disease spectra: non-DNMT3A CHIP was associated with broad multisystem disease burden, mLOY with a more focused respiratory, musculoskeletal and cardiovascular profile, whereas mLOX lacked broad age-related disease associations. Clone burden mapped to distinct proteomic programs: mLOY to neutrophil degranulation and extracellular-matrix remodeling, non-DNMT3A CHIP to myeloid immune regulation, and mLOX unexpectedly to cytotoxic lymphocyte/NK-cell responses. Mendelian randomization supported selected protein-disease relationships. Thus, age-related hematopoietic clones are not interchangeable markers of ageing but define alteration-specific systemic programs associated with distinct disease vulnerabilities.

Machine-learning models of the human microbiome are trained mostly on stool samples from single cohorts, limiting cross-body-site representation and cross-study generalization. Progress is constrained less by algorithms than by the absence of a harmonized multi-body-site corpus carrying the technical metadata needed to model, rather than ignore, batch structure. Here we release Corpusome, a harmonized two-tier cross-body-site human microbiome corpus for representation learning: a harmonized corpus of 187,546 human microbiome samples integrating standardized profiles from curatedMetagenomicData, the American Gut Project, and the EBI MGnify platform. Corpusome follows a two-tier design preserving both functional depth and cross-body-site breadth: a shotgun tier (22,588 samples, 93 studies) with species- and pathway-level profiles, and a 16S tier (164,958 samples, from a full pull of 708 MGnify studies) with genus-level profiles extending coverage to oral, skin, respiratory, and urogenital sites. It spans six body sites and two modalities, with harmonized metadata for batch-aware modelling. Body-site signal exceeds technical/source variance in the 16S tier by approximately 2.4-fold.

Single-cell technologies enable linking disease-risk variants to gene regulatory effects in specific cell-state contexts. However, most so called "single-cell eQTL" studies use a "pseudobulking" strategy to identify expression Quantitative Trait Loci (eQTLs), obscuring subtle dynamic regulatory effects of disease alleles. Here, we propose Dynema (Dynamic eQTL mapping in single cells) for fast and accurate genome-wide mapping of context-dependent and independent eQTL effects at true single-cell resolution. To identify eQTLs, Dynema uses a Poisson model with cluster robust variance estimators (CRVEs) to account for correlation of single-cell profiles from the same individual. In contrast to other common methods, Dynema achieves statistical calibration and scales to genome-wide analysis in large single-cell datasets in realistic timeframes. We applied Dynema to two independent T cell datasets and identified reproducible cell-state-dependent eQTL effects. Some cell-state-dependent eQTLs are missed by pseudobulking approaches, and many others are conditionally independent from lead eQTL effects. We show that TSPAN32 and other autoimmune loci colocalize with cell-state-dependent eQTLs. Mapping context-dependent eQTLs at single-cell resolution enables the definition of the molecular effects of complex disease alleles.

Animals rely on color to navigate complex environments, yet how eyeless organisms use chromatic information to guide food choice remains poorly understood. Here, we show that Caenorhabditis elegans exhibits robust color dependent foraging driven by microbial chromophores, preferentially consuming red while avoiding blue chromoprotein expressing bacteria across bacterial backgrounds and wild isolates. This discrimination persists in darkness and independently of photoreceptor, revealing a mechanism beyond canonical photoreception. Purified chromoproteins and bacterial metabolite fractions independently reproduce preference, demonstrating complementary chromatic and post ingestive metabolic cues. Mechanistically, blue chromoproteins generate singlet oxygen, producing oxidative stress and remodeling bacterial tryptophan and pterin metabolism, whereas red food promotes serotonin production and feeding-associated neuropeptide signaling. Disrupting serotonin biosynthesis or neuropeptide processing abolishes color preference. Together, our findings reveal a previously unrecognized, novel sensory strategy in which wavelength-selective pigment photochemistry transforms microbial color into metabolic information that is integrated through gut brain neuroendocrine signaling to guide foraging behavior in an eyeless animal.

Phage D29 infects Mycobacterium smegmatis mc2 155 and has a non-canonical lysis cassette that encodes two endolysin proteins (Lysin A and Lysin B) and a single two transmembrane domain (TMD) protein, LysA2a similar to F1 cluster phage LysF1a. A 1TMD LysF1b homolog, LysA2b, is encoded by a gene found downstream of the tape measure. Exogenous expression of both LysA2 proteins in tandem is a cytotoxic to M. smegmatis. Deletion of lysA2a produces phages that are lysis competent with a 10-minute triggering delay and 30% plaque size reduction. Deletion of lysA2b results in severe lysis defects manifest by 70% reduced plaque size, delayed lysis timing and reduced burst size. Deletion of both lysA2 genes results in phages that are viable and show lysis phenotypes like the lysF1b deletion. Genetic complementation of lysA2b deleted phage with the lysF1b gene fully complements the lysis phenotypes but alters the triggering time to that of an F1 cluster phage. Energy poisons trigger lysis prematurely in all phages with lysA2 gene deletions. Lysis recovery mutants (LRM) isolated from phages lacking the lysA2b genes generate wild type plaque size and have point mutations that map to TMD1 or the C-terminal region of the lysA2a gene. LRMs isolated from phages lacking both lysA2 genes show premature lysis and have mutations that all map to residue C31 of a novel lipoprotein (gene 64). Deletion of gene 64 does not change wild type D29 lysis phenotypes or rescue the lysis defects of any of the lysA2 mutants. A fitness/competition assay shows that loss of the lysA2 genes imposes a substantial competitive fitness cost. These finding support a lysis regulatory network model where the 2TMD protein is maintained in an inactive state until activated by its cognate 1TMD lysis regulator and the lipoprotein has accessory function that may enhance lysis efficiency.

To profile unknown ncRNAs-"dark matter" in single cells, we developed dropTotal, a high-throughput droplet-based total RNA-seq method that uses dU-modified GAT primer with temperature-ramp hybridization and droplet merge barcoding to co-detect coding and non-coding transcripts with record sensitivity (>13,500 genes/cell, including >2,000 lncRNAs and >500 sncRNAs), compatible with fresh, frozen, fixed, and FFPE tissues. Applied to ~75,000 human glioma nuclei, it captured 60,313 genes (18,681 lncRNA, 19,859 mRNAs and 6,753 sncRNAs), enabling ncRNA-driven regulatory landscape construction. In oligodendroglioma, module analysis identified recurrence-associated ncRNA-centered modules linked to therapy resistance and invasion; in glioblastoma, six cellular states showed hundreds of state-specific unannotated ncRNAs with divergent functions, from MIR222HG-mediated immune modulation to SCIRT-driven hypoxia adaptation. Alternative splicing analysis identified 428 state-specific junction markers and mapped cell-state-specific alternative splicing regulation. dropTotal offers broad application for decoding the underlying ncRNA biology and single-cell whole transcriptome regulatory mechanisms in cellular identity and disease progression.

Autosomal dominant gain-of-function mutations in platelet-derived growth factor receptor beta (PDGFRb) cause overgrowth of the skeleton and other connective tissue in Kosaki overgrowth syndrome. However, the target cell type and signaling pathways underlying PDGFRb-driven overgrowth are unknown. Normal postnatal growth is controlled by pituitary-secreted growth hormone (GH), which activates the STAT5 transcriptional factor to upregulate insulin-like growth factor 1 (IGF1). To investigate the role of the GH-STAT5-IGF1 pathway in PDGFRb-related overgrowth, we generated mice with a PDGFRb gain-of-function mutation in skeletal and fibroblast lineages, which resulted in STAT5 activation and gigantism. Conditional deletion of Stat5ab in connective tissue lineages rescued skeletal overgrowth and keloid-like fibrosis in the skin. Conditional deletion of GH receptor (Ghr) did not rescue overgrowth, indicating the physiological activator of STAT5 is not required for overgrowth. However, deletion of Igf1, the STAT5 target gene, and its receptor, Igf1r, in connective tissue, rescued the overgrowth phenotype. These findings demonstrate a GHR-independent STAT5-IGF1 signaling pathway in mutant connective tissue cells, which mediates PDGFRb-driven overgrowth in mice and potentially in humans with similar PDGFRB mutations.

Accurate determination of kinase inhibitor binding modes could provide essential information for understanding resistance mechanisms and accelerating drug discovery. While conventional structural methods such as X-ray crystallography, cryo-EM and NMR provide high-resolution information but are low-throughput and capture largely static snapshots of dynamic protein-ligand interactions Here, we introduce a single-molecule nanopore tweezer platform that functionally subtypes ATP-competitive Abl kinase inhibitors by resolving distinct ionic current signatures of Abl-inhibitor complexes. This approach distinguishes Type I, Type IIA, and Type IIB inhibitors without structural determination. We further show how clinically relevant Abl variants (T315I and E255V) reshape inhibitor engagement and binding modes. By combining baseline probability features with wavelet-based time-frequency descriptors, ensemble machine-learning models achieved 97.5% classification accuracy across seven kinase inhibitor binding modes at sub-angstrom resolution and enabled deconvolution of mixed-inhibitor samples at nanomolar concentrations. These results establish nanopore tweezers as a label-free, super-resolution platform for profiling kinase conformational states and inhibitor binding modes, complementing structural approaches and supporting precision oncology.